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ABSTRACT 

Spectroscopic analyses of Type la supernovae have shown there exist four spec- 
troscopic groups — cools, broad line, shallow silicon, and core normal — defined by the 
widths of the Si II features A5972 and A6355. 1991bg-likes are classified as "cools". 
Cools are dim, undergo a rapid decline in luminosity, and produce significantly less ^^Ni 
than normal Type la supernovae. They also have an unusually deep and wide trough in 
their spectra around 4200 A and a relatively strong Si II absorption attributed to A5972. 
We examine the spectra of supernova (SN) 1991bg and the cools SN 1997cn, SN 1999by, 
and SN 2005bl using the highly parameterized synthetic spectrum code SYNOW, and 
find general agreement with similar spectroscopic studies. Our analysis reveals that this 
group of supernovae is fairly homogeneous, with many of the blue spectral features well 
fit by Fe II. The nature of the spectroscopic commonalities and the variations in the 
class are discussed. Finally, we examine intermediates such as SN 2004eo and discuss 
the spectroscopic subgroup distribution of Type la supernovae. 

Subject headings: supernovae: individual (1991bg, 1997cn, 1999by, 2005bl) 



1. Introduction 

Type la supernovae (hereafter SNe la) are important objects in the study of nucleosynthesis, 
stellar evolution, and modern cosmology. It is believed that SNe la are thermonuclear explosions 
of carbon-oxygen white dwarfs which have approached to within ab out 1% of the Chandrasekhar 



Mass (1.39Mq) (iHillebraiidt Niemevedl2000l : iHoflich et al.ll2010l ). In cosmology, SNe la were 
suggested early on ( Baadel 1938 ) to be useful as standard candles for making accurate distance 



measurements in the determination of cosmological parameters, because of both their significant 
homogeneity in absolute bolometric magnitude and large apparent brightness. With modern su- 
pernova surveys such as the Supernova Cosmology Project, ESSENCE Supernova Survey, Nearby 
Supernova Factory, Palomar Transient Factory, Supernova Legacy Survey, and the Large Synoptic 
Survey Telescope, accurate measurements of these parameters based on SNe la at varying redshifts 
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Law et alj l2009l : IConlev et al.l 1201 ll : ILSST Science Collaborations et al 



have been and c ontinue to be mad e (IPerlmutter et alJll999l: iKrisciunas l2008l : ICopin et alJ 120061 : 



200 



It has even been 

possible to measure the rate of expansion of the universe and put limits on cosmoloKical models 



(iRiess et al 



1998 



Goobar et al 



200C; 



Wood-Vasev et alJl2007l : iRiess et alJl2009l : Ifficken et alJl2009l : 



Chuang &: Waneil201ll : [Parkinson et al.ll20ld ). It is because of this important role SNe la play that 
a deep understanding of their properties, variations, and evolution is vital. 

In our study we have narrowed our focus to the "cool" spectroscopic subgroup, epitomized 
by SN 1991bg which occurred in NGC 4374, and so-called "intermediates" which could potentially 
bridge the gap between the four spectroscopic subgroups. We used the highly-parameterized syn- 
thetic spectrum code SYNOW t o model the spectra of several co ol SNe la. To better understand 



the SYNOW code and its use see Branch et al. (2007 



2003 



2002) , and for more in depth technical 



details see iFished ()2000l ). The spectra range from a few days before B maximum to as late as one 
month post B maximum in the optical spectra with wavelengths from the Ca II H&K feature in the 
blue to the Ca II infrared triplet. The spectra have been corrected for redshift of their host galaxies 
and normalized to remove the slope of the continuum, as we are prima rily interested in sp ectral 
features for making line identifications. This normalization, described in I Jefferv et al.l (j2007l ). also 
serves to facilitate the comparisons between different supernovae. Extremely noisy spectra have 
been smoothed using boxcar smoothing. 

In f|2]the method of spectroscopic sub-classification is discussed. In f|3]the properties, progeni- 
tors, and possible explosion mechanisms of 1991bg-likes are reviewed, and in ^our line identifica- 
tions are given and compared with similar studies. In fj5]we discuss the potential of intermediates 
to fill in the gaps between the seemingly discrete subgroups. 



2. Spectroscopic Subgroups 



For decades, SNe la had been accepted for use as standard candle s based on their 



ent homogeneity, but i n 1991 two very peculiar supernovae, SN 1991bg (IFilippenko et al. 



Leibundgut et al.l Il993l ) and SN 1991T (IWaagen et al.l Il99ll: IFilippenko et al.l Il992bl ) challenged 



the validity of such use. According to work completed by iLi et al 



2001 



appar- 



1992a 



and recently updated 

(|Li et al.ll2010l ) the peculiarity rate for SNe la could be as high as 30%, with rates of 1991bg-likes 
up to 15%. These observed peculiarities instigated the search for parameters which could be used 
to find a photometric and spectroscopic standard relation, in order to make SNe la standardizable 
candles. 



The first relation was found by iPhillipsI (j 19931 ) who used a small sample of well-observed 
SNe la to identify a very strong correlation between peak luminosity and initial decline rates, 
parameterized by Amis. The Amis parameter measures the the decline in brightness of the B 
band from maximum until 15 days after maximum. The Phillips relation indicates brighter objects 
will have slower decline rates than dim objects which decline rapidly. An analogous spectroscopic 
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relation was suggested by iNugent et alj (jl995l ). Nugent's spectral relation measures the ratio of 
the strength of the Si II features A5972 and A6355, denoted by 7^(Si II). The standard relations, 
both photometric and spectroscopic, are due in large part to temperature differences caused by the 
amount of radioactive ^^Ni synthesized during the explosion. 

Single parameter modeling of SNe la is indeed useful, but in cases of extremely peculiar SNe, 
such as SN 1991bg and SN 1991T, other means of accurate and efficient sub-classificatiqn are 
necessary. Several such classification schemes exist. One scheme proposed bv lBenetti et al.l (120051) 



is bas ed on the evolution of Si II A6355 line velocity, while another, proposed by [Branch et al 
(|2006l ) uses the pseudo-equivalent widths of the Si II lines A5972 and A6355. 



By measuring the velocity of Si II A63 55 10 days past B max as well as the evolution of the 
velocity of this feature iBenetti et al.l (|2005l ) found not only a spread in velocities between different 
SNe la, but also a distribution i n velo city gradients. Plotting these gradients reveals three groups 
of SNe la which iBenetti et al.l (|2005l ) labels as f aint (FAINT) , high temporal vel ocity gradient 
(HVG), and low temporal velocity gradient (LVG) (jBenettilboOsI : benetti et ahlboosl ). The FAINT 
subgroup corresponds to subluminous 1991bg-like events with low expansion velocities and a large 
velocity gradient. The HVG group contains normal SNe la with high expansion velocities and a 
high velocity gradient. Finally, the LVG group contains both normal SNe la as well as the brightest 
SNe la. This group has, on average, lower expansion velocities than the HVG group and a low 
velocity gradient. 



In the scheme of [Branch et al.l (|2006l ). the pseudo-equivalent widths of Si II lines A5972 and 
A6355 are measured and plotted (since emission lines are present in the spectra, these are not 
true equivalent widths but rather widths whose limits of integration have been chosen by eye). 
We followed the example of [Branch et al.l (|2006l ) and plotted the width of the A5972 feature vs. 
the width of the A6355 feature for a sample of SNe la. This plot can be seen in Figure [T] and 
shows a cluster of core-normal (CN) SNe la with three discrete branches representing the peculiar 
spectroscopic subgroups of cool (CL), broad-line (BL), and shallow silicon (SS). The CLs correspond 
to Benetti's FAINT, BL tend to fall into the HVG group, and the LVG contains both CN and SS. 

In this method, each of the spectroscopic subgroups have their own differing properties and 
are named accordingly. CNs are the most common and well studied SNe la and are found in the 
region of highest density in the pseudo-equivalent width (WW) plot of Figure [H SN 1994D is an 
example of a CN SNe la. BL SNe la have Si II lines whose absorptions are broader and deeper 
than CNs and thus form the lower right branch of the WW plot in Figure [H The SNe la in the 
BL subgroup generally have th e same ions that ar e present in CN, but with higher velocities. One 
such BL SNe la is SN 2002bo (|Branch et ahlboOfil ). The SS SNe fall to the lower left of the plot in 
Figure [1] because of their smaller values for Si II widths, but otherwise they have normal SNe la 
spectra. (The CLs will be discussed in depth in ^) Although the plot indicates there are separate 
spectroscopic subgroups the borders are not defined by a clear demarcation; instead some overlap 
exists. We call the SNe la which lie between different subgroups "intermediates". 
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3. Cool Type la Supernovae 



Spectroscopically peculiar SNe la which fall into the cool subgroup share several distinct charac- 
teristics. Most notably, cools are subluminous. SN 1991bg had V and B maxima wh ich were 1.6 and 
2.5 magnitudes, respectively, lower than ordinary SNe la (jFilippenko et al.lll992al ). The CLs also 
have much faster declines from peak luminosity than do CNs. Typically, CLs have Amir^{B) ~ 1.9, 
while CNs have Ami5{B) 1.1 (ITaubenberger et al.l l2008l ). Light curves of SN 1991bg and 
SN 1994D can be found in Figure [2] showing appropriate Amis (5) values. It has been estab- 



lished that ^^Ni production is fu ndamental to determining SNe la peak luminosities ( Arnett et al 



19851 : iBranch &: Tammannlll992l ). Analyses of SN 1991bg and other CLs ha ve indicated ^^Ni pro - 
duction of only about 1/6 the amount synthesized during a normal event (jMazzali et al.l 119971 ). 
Spectroscopically, we see primarily that CLs have an unusually deep and wide absorption trough at 
around 4200 A and a very strong Si II A5972 absorption as shown in Figure [3 The 4200 A trough 
is well fit by Ti II which can be explained by its low excita tion temperature and the cool nature 
of these SNe la ([Filippenko et al.lll992al : iMazzali et al.lll997l ). Finally, CL s also have a strong and 
narrow absorption at about 5675 A usually attributed to Na I D A5893 ( Filippenko et al.lll992al ; 
Garnavich et al.ll2004l ). 



3.1. Progenitors 



To date, many theories and models have been proposed in attempts to describe the progenitors 
which lead to SNe la, as well as to account for their diversity. Two common scenarios leading 
to SNe la are genera lly accepted. One scenario is a single degenerate (SD) case proposed by 
Whelan IbenI (jl973l ) in which a white dwarf (WD) accretes mass until the Chandrasekhar Mass 
{Mch) is approached, igniting carbon in the core, undergoing thermonuclea r instability, and ending 
with the c omplete d i srupt ion of the WD. A competing scenario proposed bv llben &: Tutukovl (jl984l ) 
as well as IWebbinkI (| 19841 ) involves two orbiting WDs, where loss of angular momentum due to the 
emission of gravitational waves leads to an inspiral and eventual merger. The merger of two WDs 
is referred to as double degenerate (DD). Extensive work has been done to investigate how either 
a SD scenario or a DD scenario can possibly explain both the homogeneity of SNe la while at the 
same time account for the observed diversity. 

SNe la appear to be homogen eous because the str ucture of the WD progenitor and the explosion 
is determined by nuclear physics (iHofiich et al.ll2010l ). That is, the WD is supported by degenerate 
electron pressure, thermonuclear energy gives the energy production during explosion^ and light 
curves are driven by the decay of synthesized radioactive ^^Ni. According to [Hoflichj ( 20061 ). the 
aspects most important to understanding the diversity are the conditions just prior to explosion and 
the type of burning mechanism. The basic types of burning can be distinguished as detonation and 
deflagration. Deflagration occurs with a burning fr ont propagating due to hea t transport across the 
front at a speed of 0.5 — 0.8% the speed of sound (jTimmes &: Wooslevlll992l ). On the other hand. 



- 5 - 



in detonations compressiona l heat ignites carbon and oxygen in front of a shock which propagates 
supersonicaUy (jHofiichl 120061 ) . 



SD models using delayed detonation — a process where deflagr ation occurs follo wed by a tran- 
sition to detonation — show promise in reproducing observations (iKhokhlovi Il99ll ) . Varying the 
amount of pre-expansion and burning during the deflagration phase as well as a few other sec- 
ondary parameters can expl a,in the diversi t y am ong SNe la in a SD scenario (jHoeflich et al.l 119931 : 
Hofiich et al.ll20ld ) . Namely, iHoflich et al.l (j2010l ) decreases the density where transition from defla- 
gration to detonation occurs and in so doing reduces the amount of ^^Ni synthesized thus producing 
subluminous SNe la. It should be noted that the trigger for transition from deflagration to det- 
onation is still not fully understood. When considering the DD scenario recent work using the 
smoothed particle hy drodynamics code GADGET S shows promise in reproducing subluminous SNe la 



akin to 1 99 Ibg- likes (jPakmor et al 



2011 



20101 ). They find that this type of merger reproduces the 



observed results and may be a path to the production of 1991bg-like SNe la. 

Another method used to investigate SNe la progenitors is population synthesis. Through this 
method, the observed delay time distribution (DTD) of SNe la can be used to place limits upon prq - 
genitor models and formation scenarios. Taking an observational approach iMaoz &: Badened (120101) 



Maoz fc Badenes 



derive the SN rates and DTD using supernova remnants in the Magellanic Clouds. 
(2010) find that DTD is prop o rtiona l to in agreement with the DD DTD found from population 
synthesis of DD. iMaoz et al.l (|201ll ) also confirmed this DTD result using star formation history. 



The general findings of popu lation synthesis ar e that neither channel, SD or DD, can produce th e 



observed number of SNe la (Ruiter et al 



200S 



Mennekens et al.ll20ld ). iMennekens et all ()20ld ) 



using the Brussels population synthesis code, find that the morphological shape of the observed 
DTD cannot be fit by the SD scenario alone. Instead, a combination of both SD and DD scenar- 
ios will fit the morphological shape but produce s approximately three times fewer SNe la than is 
observed. A similar study bv lRuiter et al.l (j2009l ). using the population synthesis code StarTrack, 
also find that the combination of SD and DD scenarios fits the observed shape of DTD, but ac- 
counts for ten times fewer SNe la than observed. Population synthesis would then seem to indicate 
that neither scenario alone can produce the morphological shape of the DTD and that even the 
combination of both cannot produce the absolute number of observed SNe la. 



4. Line Identifications in Cools 

We begin our analysis of CL SNe la by making line identifications. Identifying lines in SNe 
can be a difficult process and is further complicated in CLs due to the large number of ions with 
lower excitation temperatures and many metals which potentially are present as blends in the 
4200 A trough. Our sample contains spectra from SN 1991bg, SN 1997cn, SN 1999by, and SN 
2005bl. We find the nine ions— O I, Na I, Mg II, Si II, S II, Ca II, Ti II, Cr II, and Fe II— of 
varying strengths and levels of stratification fit most of the observed spectral features well for all 
epochs considered with a tenth ion, Ca I, present in the early spectra. 
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All with respect to B maximum, we will make comparisons of five subsets of spectra from dif- 
ferent epochs. In S i4.1l we compare SN 1999by and SN 2005bl at day —3. The observed and synthetic 
spectra for this comparison are shown in Figure |4] and the fitting parameters are listed in Tables [TH2j 
In SJwe analyze the spectra of SN 1991bg, SN 1997cn, SN 1999by, and SN 2005bl just past B max. 
These spectra can be found in Figures [5H6] with fitting parameters listed in Tables [SHSl Synthetic 
spectra and line identifications for SN 1991bg at day 18 and SN 2005bl at day 19 are summarized in 
§4.3l and displayed in Figure [8] with fitting parameters listed in Tables [THHl In §4.41 spectra from ap- 
proximately one month post B maximum for SN 1991bg, SN 1997cn, and SN 1999by are discussed. 



9HT31 A great deal of spectroscopic analysis has been performed on SN 1991bK bv FilioDenko et al. 


( 1992a 


): Leibundsut et al. (1993): 


Turatto et al. (1996): Mazzah et al. C 


L997) and ] 


Branch et al. 


(2006, 


2008. 


2009). Turatto et al. 


(1998) have analvzed SN 1997cn. Garnavich et a] 


. (2004 


) and 


Hoflich et al. 


(2002) considered SN 1999bv. Taubenbereer et al. (2008) and 


Hachineer et al. ( 


2009) 



performed an analysis of SN 2005bl. 



4.1. Day 3 Pre B Maximum Spectra 



The day —3 spectra are presented in Figure HI The SYNOW spectra for SN 1999by and 
SN 2005bl are overplotted with their respective observed spectra, and spectral absorptions are 
labeled. The synthetic spectra parameters are given in Tables [TH2j 

The broad absorption trough at 4200 A, characteristic of the CL subgroup, is well fit by a 
combination of Ti II and Mg II. This tr ough is continually fit by both Ti II and Mg II in all 
modeled spectra. iGarnavich et al.l (120041 ). while modeling SN 1999by with SYNOW, find the 
absorption at 5000 A to be fit by Mg I, whereas we attribute this feature to Ti II and find no need 
of Mg I. Both Ti II and Mg I have an absorption at about 5000 A; however, SYNOW fits with Ti II 
and without Mg I give accurate fits for this absorption. Mg I, while possibly present, when included 
is not the main contributor to this feature and has no other absorptions or emissions which can 
clearly be attributed to it. Thus we consider the identification of Mg I to not be required. Ti II 
also contributes to the 5780 A feature. The 7630 A feature, blended into the O I triplet, belongs 
to Mg II. For these reasons Ti II and Mg II can be considered definite. Other definite ions include 
Si II, Ca II, Fe II, and O I which all contribute to multiple absorptions. 

The other four ions, Ca I, Cr II, S II, and Na I, only contribute to one absorption each and, 
although quite likely, are not d efinite. The "W" fea. ture at 5330 A, often attributed to S II, may also 
be fit by features due to Sc II; [Branch et al.l (|2006l ) have indeed produced SYNOW fits with this 
attribution. Although synthetic fits with Sc II are possible, S II is more likely in the SNe la scenario, 
we discuss this further below. Na I creates a shoulder in the 5650 A emission, and although very 
weak in the early spectra, becomes much more conspicuous later. Based on the evolution of multiple 
epochs, even though Na I is wea k early, it is highly likely to be present. Na I is a common identifica- 



tion of this feature in CL SNe la (jFilippenko et al.lll992al : IGarnavich et al.ll2004l : lTaubenberger et al 
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20081 ) : however, iMazzali et alj (|l997l ) when using W7 cannot repr oduce this feature withou t adding 



extra Na I and modifying their abundance distributions, while iLeibundgut et al.1 (|l993l ) identify 
it as a blend of Si II. Ca I matches the absorption a t 5980 A in these spec tra, but disappears 
sometime after day 4 post B max and before day 18. iGarnavich et al.1 (120041 ) identifies two Ca I 
absorptions as late as day 7 in SN 1999by and iTaubenberger et al.l (j2008l ) make no Ca I identifica- 
tions for SN 2005bl. The remaining ion, Cr II, is fo und at the 4700 A a b sorpt ion and further helps 
to block the emission a t 4600 A, in agreement with Taubenberger et al. (j2008l ). but unidentified by 



Garnavich et al.l (|200J). 



Aside from ICarnavich et al.l (|2004l ) identifying Mg I in SN 1999by at 5000 A instead of Cr II, 
our ot her line identifications are in good agreement. An analysis performed by iHachinger et al 



([20091) on SN 2 005bl identifies sma l l amo unts (less than 10% of the mass fraction) of C II in 
agreement with ITaubenberger et al.l (j2008l ) who also find C II. C II i s possi ble, however does not 
add any benefits to our synthetic fits and as noted by IHachinger et al.l (j2009l ) is a break from other 
CL SNe la such as SN 1991bg and SN 1999by in which C II has not been iden tified. In fact, C II 



is more ofteri associated with superlurn inous, rather than subluminous, SNe la (jHowell et al. 



2006; 



Hicken et al. 


2007; 


Parrent et al. 


2011 



4.2. Days 2, 3, and 4 Post B Maximum Spectra 

Figures [5H6] show the observed and synthetic spectra for SN 1991bg, SN 1997cn, SN 1999by 
and SN 2005bl at days 2, 3, 3, and 4, respectively. The spectrum for SN 1997cn ends at 7400 A, 
which is before the O I feature and Mg II absorption at 7630 A; however these ions are still used 
in the SYNOW fits as they contribute to the overall shape of the spectra. The fitting parameters 
used for these synthetic spectra are given in Tables [SHSl 

As usual we consider Si II, Ca II, Fe II, and O I definitely present for their contribution to 
multiple spectral features. Also, Ti II and Mg II are definite ions except in the case of 1997cn 
where the Mg II feature at 7630 A is not available to be modeled. Due to the attribution of Mg II 
to the absorption trough in SN 1997cn we still consider it highly likely. Again, Cr II continues to 
be present at 4700 A and helps to reduce the emission just blueward of its absorption. We consider 
the presence of Cr II not only plausible but likely. Further, during these early spectra Ca I fits the 
absorption at 5980 A. 



In comparison with the analysis by iTuratto et al.l (|l998l ) on SN 1997cn, we continue to see 



good agreement wit h the primary differences being Ca I and Cr II which they do not identify. 



Turatto et al.l ( 19981 ) also identifies Ni II at 4067 A, which is possible, but we find no benefit in the 



inclusion of Ni II in our fits. 

Close examination of the fitting parameters given for these SNe in Tables [3H6] will show that 
both 1997cn and 2005bl have photospheric velocities (vphot) of 7, 600 km s~^ which is significantly 
lower than 1991bg at 11, 000 km and 1999by at 10, 000 km s"'^. Looking forward to days 18 and 
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19 we see that 2005bl has a lower gradient in fp^ot than does 1991bg and differs by only 100 km s 
By days 28, and 29 the differences in Vphot between the SNe are insignificant. 

Figure [7] compares the fit of SN 1999by at day +3, using S II and Sc II to fit the "W" feature 
usually attributed to S II. It is hard to definitively decide between the two fits and without the 
strong S II A6355 line, the identification of SN 1999by as a SNe la becomes weaker. This effect is 
important but needs to be followed up in future work with more detailed spectral modeling. 

4.3. Days 18 and 19 Post B Maximum Spectra 

All the same ions as discussed in §SilHll2l except Ca I, are still present by days 18 and 19, 
but the spectra have several noticeable differences from the early time spectra. (Synthetic and 
observed spectra for SN 1991bg and SN 2005bl can be seen in Figure [8j) Primarily the differences 
include the ions having deeper absorptions, higher emissions, and lower Vphot- The most obvious 
difference is the Na I absorption at 5700 A, which is no longer just a shoulder in an emission but a 
rather strong and sharp absorption. We still find Si II, Ca II, Fe II, Ti II, Mg II, and O I definite 
whereas Na I is now much more likely. Cr II and S II are plausible for the same reasons as before 
and Ca I is no longer used in the synthetic spectra. 

4.4. Days 28, 29, 31, and 32 Post B Maximum Spectra 

Figures [9HT0] have the synthetic and observed spectra for SN 1991bg, SN 1997cn, and SN 1999by 
where the absorptions are well fit; however, in the red the synthetic emissions are weak. This lack 
of strong emissions can be explained by the assumptions of SYNOW itself. SYNOW assumes a 
perfect blackbody emitter and absorber with a sharp photosphere surrounded by a homologous and 
spherically expanding ion cloud where lines form by resonance scattering treated in the Sobolev 
approximation. These assumptions do not allow for net emission profiles. Once ions begin to enter 
net emission SYNOW can no longer reproduce the strong emission peaks well. To overcome this 
issue it is possible to increase greatly the optical depth and rescale, but better fits are often obtained 
simply by leaving these emissions alone and focusing on the absorptions. 

The same nine ions are present in these spectra as the previous, but O I has been strongly 
stratified to fit the three separate absorptions at about 7530 A. Furthermore, a deep and wide 
absorption around 6700 A to 7000 A has appeared. Despite our best efforts we were unable to fit 
this trough in a reasonable way. We tried a variety of ions including C I and O II, which could 
fill some of the trough but not without heavy restriction in Vmin and Vmax ■ None of the attempts 
proved successful and we do not find any of these ions likely. 
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5. Intermediate Type la Supernovae 



Figure [TT] is another WW plot very similar to Figure [H but in tliis WW plot, intermediates are 
included. Intermediates are SNe la that do not necessarily fall into any particular subgroup; instead 
they fall in between, sharing properties of the different subgroups. What we find when including the 
intermediates is that the subgroups begin to blend, indicating that rather than being discrete, SNe la 
subgroups, as delineated by Branch, are a continuous distribution based on multiple parameters. 
Branch et al.l (|2006l . 120091 ) did extensive work on spectroscopic subclassification of SNe la and came 
to a similar conclusion. 



Recent work of iMaeda et al.l (|20ld ) indicates that some of the spectral diversi ty of SNe la 



i s a r esult of the viewing angle on an asymmetrically exploding WD progenitor. iMaeda et al 
(j2010l ) used the Benetti classification scheme of FAINT, HVG, and LVG. Using the emission lines 
Fe II A7155 and Ni II A7378 lines, which occur during the nebular phase, they measured the Doppler 
shift in the ejecta. Blue-shifted lines will be on the near side and red-shifted will be on the far side 
of an exploding SN. They found a diversity in the velocities of the blue-shifted and red-shifted lines 
implying that the initial SN ignition occurred off-center. Measurement of red-shift and blue-shift 
in HVG and LVG SNe la show that HVGs are preferentially red-shifted and LVGs blue-shifted. 
Maeda et al.l (|2010l ) conclude that LVGs are viewed in the direction of the initial spark and HVGs 
are viewed opposite the spark. Statistical treatment of this result indicates a high pr obability that 



the di fferences of the two groups is entirely an effect of viewing angle. The work of iMaeda et al 



(|2010l ) explained why the differences between the HVG and LVG groups occurred, but did little to 
address the FAINT group which they acknowledge may arise from a completely different explosion 
mechanism then HVGs and LVGs. 



6. Discussion 



Due to the high level of importance placed on SNe la, a complete understanding of their 
physical nature is paramount. One method to probe this nature is examination of their differences. 
Thus, we analyzed the spectroscopically CL subgroup of SNe la. In summary, SNe in this subgroup 
are underluminous, rapid decliners, and poor ^^Ni producers with deep and wide absorption troughs 
around 4200 A due to Ti II and unusually strong Si II A5972 absorptions. The path leading to 
subluminous SNe la is presently uncertain. SD scenarios in which a delayed-detonation, where 
a subsonic deflagratio n front transitions into a detonation, are strong theoretically and account 
well for obse rvations (|Hofiich et al.ll2010l ). There also exist competing DD scenarios such as that 



presented by lPakmor et al 



(2011 



which indicate SNe la potentially originate from the merger of 
two CO WDs. Each scenario has its own strengths and weaknesses and the question of progenitors 
remains open. Furthermore, population synthesis studies tell us that neither scenario alone can 
account for the total number of observed SNe la. 



We have modeled the spectra of several CL SNe la with the fast and highly parameterized 
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synthetic spectrum code SYNOW and compared our findings with those of similar studies. We 
find good agreement with other studies. Our work indicates many ions are potentially present in 
CLs, and line identifications are complicated for this reason. We show that Si II, Ca II, Fe II, Ti II, 
Mg II and O I are the most definite ions present in the spectra of cools. It is also very likely that 
Na I is present. Three other ions, Ca I, Cr II, and S II, are considered as possibilities. S II is quite 
likely, but Ca I and Cr II only occur in one location in our synthetic spectra and they have not been 
definitively explained in other studies. Therefore we consider these ions possible but not definite. 

Examining the velocity evolution of Si II A6355 shows the primary difference among this 
subgroup. From our models we see that SN 2005bl and SN 1999by have lower Si II A6355 velocities 
than SN 1991bg and SN 1997cn in early time spectra but have lower velocity gradients so that by 
about day 18 post B maximum all four CLs have approximately the same velocity. Si II A6355 
velocities are plotted in Figure [12] for several CLs and CNs. Comparing our velocities with those 
of other studies show insignificant differences. Overall, the appearance of the velocity distribution 
for all SNe la is quite continuous. 

Our examination of this subgroup indicates that 1991bg-likes are fairly homogeneous amongst 
themselves (albeit somewhat due to construction) and not only share the same ions, but fitting 
parameters and properties as well. Upon consideration of intermediates the distinctions among the 
subgroups begin to blur even more and the distribution of SNe la seems increasingly continuous. 
With a large and well-observed sample this distribution may take a form similar to the hypothetical 
distribution displayed in Figure [HI 

A continuous distribution of SNe la could have regions of high density and low density in 
the WW plots, but the difference in density should be fairly small. These dense regions would 
contain the subgroups. The highest density region would contain the most common SNe la, the 
CNs. The lowest density regions would contain SNe la which are steps away from one subgroup 
toward another; intermediates. The reason for different regions of high density and low density, 
i.e. the subgroups, could indicate sets of secondary parameters which are slightly more stable and 
thus more likely to be reproduced for WDs transitioning into SNe la, rather than fundamentally 
different processes. 

Present work on the diversity of SNe la is incomplete but is leading to the conclusion that 
rather than being made of discrete subgroups, SNe la are a continuous distribution of the same 
phenomenon. However, for more definitive conclusions it is necessary to have a much larger and 
well-observed sample coupled with more robust explosion models. 
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Fig. 1. — Pseudo-Equivalent Widths of A5972 vs. A6355. Data comes from iBranch et al.l ( 2003 ) as 
well as private communication with D. Branch. 
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Fig. 2. — CN vs. CL Decline Rate Comparison s. Ami 5 valu e s orig inate from iPastorello et al. 
(j2007l ) , light curv e data for 1991bK comes from iTuratto et al.l (|l996l ) , and light curve data for 
1994D comes from lPatat et al.l (|l99d ). 
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Fig. 3. — CN vs. CL Spectra Comparisons 
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Fig. 4. — 2005bl and 1999by Day 3 Pre B Maximum. The telluric oxygen feature at 7594 A is not 
marked. 
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Table 1. Fitting Parameters for 1999by Day 3 Pre B Maximum. 
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*A11 velocities are given in 1000 km s ^ 



Table 2. Fitting Parameters for 2005bl Day 3 Pre B Maximum. 
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*A11 velocities are given in 1000 km s 
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Fig. 5. — 1997cn Day 3 Post B Maximum and 1991bg Day 2 Post B Maximum. The telluric oxygen 
feature at 7594 A is not marked. 
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Table 3. Fitting Parameters for 1991bg Day 2 Post B Maximum. 
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*A11 velocities are given in 1000 km s 



Table 4. Fitting Parameters for 1997cn Day 3 Post B Maximum. 
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*A11 velocities are given in 1000 km s 
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Fig. 6. — 1999by Day 3 Post B Maximum and 2005bl Day 4 Post B Maximum. The telluric oxygen 
feature at 7594 A is not marked. 
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Fig. 7. — 1999by Day 3 Post B Maximum with the feature usually attributed to S II fit with both 
S II and Sc II. The telluric oxygen feature at 7594 A is not marked. 
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Table 5. Fitting Parameters for 1999by Day 3 Post B Maximum. 
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40 


O I 














Na I 


0.5 


7000 


10 


1.0 


12.5 


13.5 


Mg II 


40.0 


7000 


10 


1.0 





40 


Si II 


100 


7000 


10 


0.5 





40 


S II 


5.0 


7000 


10 


0.2 





40 


Ca I 


12.0 


7000 


10 


1.0 





40 


Ca II 


1100 


7000 


10 


1.3 





40 


Ti II 


6.0 


7000 


10 


2.0 





40 


Cr II 


1500 


7000 


10 


1.0 





40 


Fe II 


4.0 


7000 


10 


1.0 





40 



*A11 velocities are given in 1000 km s 



Table 6. Fitting Parameters for 2005bl Day 4 Post B Maximum. 



Element 


T 


Tea:c (K) 


Vphot 


Ve 






HVe O I 














I 


4.5 


7000 


7.6 


1.0 





40 


Na I 


0.1 


7000 


7.6 


1.0 


10 


12 


Mg II 


10.0 


7000 


7.6 


1.0 





40 


Si II 


10.0 


7000 


7.6 


0.6 





40 


S II 


1.0 


7000 


7.6 


0.2 





40 


Ca I 


8.0 


7000 


10 


1.0 





40 


Ca II 


800 


7000 


7.6 


1.3 





40 


Ti II 


7.0 


7000 


7.6 


2.0 





40 


Cr II 


600 


7000 


7.6 


1.0 





40 


Fe II 


6.0 


7000 


7.6 


1.0 





40 



*A11 velocities are given in 1000 km s 
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Fig. 8.— 2005bl Day 19 Post B Maximum, and 1991bg Day 18 Post B Maximum. 
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Table 7. Fitting Parameters for 1991bg Day 18 Post B Maximum. 



Element 


T 


Te^c (K) 


Ujjhot 








UVe O I 


2.0 


7000 


8 


1.6 





40 


O I 


0.7 


7000 


8 


1.0 


15.2 


16 


Na I 


9.0 


7000 


8 


0.5 


11 


12 


Mg II 


3.5 


7000 


8 


1.0 


12 


40 


Si 11 


50 


7000 


8 


0.5 





40 


S II 


0.8 


7000 


8 


1.0 


14.5 


40 


Ca II 


4000 


7000 


8 


1.3 





40 


Ti II 


12 


7000 


8 


1.0 





40 


Cr II 


2500 


7000 


8 


1.0 





40 


Fe II 


70 


7000 


8 


0.7 





40 



*A11 velocities are given in 1000 km s 



Table 8. Fitting Parameters for 2005bl Day 19 Post B Maximum. 



Element 


T 


Te^c (K) 


fphot 








HVi O I 


1.0 


7000 


7.5 


1.6 





40 


O I 














Na I 


3.0 


7000 


7.5 


0.5 


11 


12 


Mg II 


2.0 


7000 


7.5 


1.0 


12 


40 


Si II 


40 


7000 


7.5 


0.5 





40 


S II 


0.7 


7000 


7.5 


1.0 


14.5 


40 


Ca II 


4000 


7000 


7.5 


1.0 





40 


Ti II 


17 


7000 


7.5 


1.0 





40 


Cr II 


2800 


7000 


7.5 


1.0 





40 


Fe II 


30 


7000 


7.5 


0.8 





40 



*A11 velocities are given in 1000 km s 
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Fig. 9.— 1997cn Day 28 Post B Maximum, 1999by Day 29 Post B Maximum, and 1991bg Day 29 
Post B Maximum. 
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Table 9. Fitting Parameters for 1997cn Day 28 Post B Maximum. 



Element 


r 


Tea:c (K) 


Vphot 


Ve 




Vmax 


HVe O I 


0.6 


7000 


7.4 


1.6 





40 


O I 


0.4 


7000 


7.4 


1.0 


14 


16 


Na I 


6.0 


7000 


7.4 


0.7 


11 


40 


Mg II 


7.0 


7000 


7.4 


1.0 


12 


13 


Si II 


10 


7000 


7.4 


0.5 





40 


S II 


0.4 


7000 


7.4 


1.0 


14 


40 


Ca II 


4000 


7000 


7.4 


1.3 





40 


Ti II 


8.0 


7000 


7.4 


1.0 





40 


Cr II 


3000 


7000 


7.4 


1.0 





40 


Fe II 


50 


7000 


7.4 


0.6 





40 



*A11 velocities are given in 1000 km s ^ 



Table 10. Fitting Parameters for 1991bg Day 29 Post B Maximum. 



Element 


T 


Texc (K) 


Vphot 


Ve 




Vmax 


UVe O I 


0.6 


7000 


7.6 


1.6 





40 


O I 


0.4 


7000 


7.6 


1.0 


14 


16 


Na I 


6.0 


7000 


7.6 


0.7 


11 


40 


Mg II 


7.0 


7000 


7.6 


1.0 


12 


13 


Si II 


10 


7000 


7.6 


0.5 





40 


S II 


0.4 


7000 


7.6 


1.0 


14 


40 


Ca II 


4000 


7000 


7.6 


1.3 





40 


Ti II 


8.0 


7000 


7.6 


1.0 





40 


Cr II 


3000 


7000 


7.6 


1.0 





40 


Fe II 


50 


7000 


7.6 


0.6 





40 



*A11 velocities are given in 1000 km s 



Table 11. Fitting Parameters for 1999by Day 29 Post B Maximum. 



Element 


T 


Tea:c (K) 


Vphot 


Ve 


Vmin 


Vmax 


HVe O I 


0.6 


7000 


7.6 


1.6 





40 


O I 


0.4 


7000 


7.6 


1.0 


14 


16 


Na I 


3.0 


7000 


7.6 


0.7 


11 


40 


Mg II 


7.0 


7000 


7.6 


1.0 


12 


13 


Si II 


15 


7000 


7.6 


0.5 





40 


S II 


0.4 


7000 


7.6 


1.0 


14 


40 


Ca II 


4000 


7000 


7.6 


1.3 





40 


Ti 11 


8.0 


7000 


7.6 


1.0 





40 


Cr II 


3000 


7000 


7.6 


1.0 





40 


Fe II 


50 


7000 


7.6 


0.6 





40 



*A11 velocities are given in 1000 km s 
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Fig. 10.— 1999by Day 31 Post B Maximum and 1991bg Day 32 Pre B Maximum. 
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Table 12. Fitting Parameters for 1999by Day 31 Post B Maximum. 



Element 


T 


Te^c (K) 


Uphot 








UVe O I 


0.6 


7000 


7.5 


1.6 





40 


O I 


0.1 


7000 


7.5 


1.0 





40 


Na I 


3.5 


7000 


7.5 


0.6 


12 


40 


Mg II 


0.1 


7000 


7.5 


1.0 


12 


13 


Si 11 


10 


7000 


7.5 


0.6 





40 


S II 


0.4 


7000 


7.5 


1.0 


14.5 


40 


Ca II 


4000 


7000 


7.5 


1.3 





40 


Ti II 


7.0 


7000 


7.5 


1.0 





40 


Cr II 


3300 


7000 


7.5 


1.0 





40 


Fe II 


35 


7000 


7.5 


0.8 





40 



*A11 velocities are given in 1000 km s 



Table 13. Fitting Parameters for 1991bg Day 32 Post B Maximum. 



Element 


T 


Te^c (K) 


^phot 








HVi O I 


0.6 


7000 


7.6 


1.6 





40 


O I 


0.5 


7000 


7.6 


1.0 


14 


16 


Na I 


6.0 


7000 


7.6 


0.6 


11 


40 


Mg II 


4.0 


7000 


7.6 


1.0 


12 


14 


Si II 


10 


7000 


7.6 


0.5 





40 


S II 


0.1 


7000 


7.6 


1.0 


14.5 


40 


Ca II 


4000 


7000 


7.6 


1.3 





40 


Ti II 


7.0 


7000 


7.6 


1.0 





40 


Cr II 


3500 


7000 


7.6 


1.0 





40 


Fe II 


40 


7000 


7.6 


0.6 





40 



*A11 velocities are given in 1000 km s 
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Fig. 11. — Pseudo - Equiv alent Widths of A5972 vs. A6355 including intermediates. Data comes 
from [Branch et al.l (|2009l ) as well as private communication with D. Branch. 
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Fig. 12. — Si II A6355 velocities. Filled symbols are CL SNe la and diamond s ymbo ls represent 
veloci ti es from our data . Oth er li ne velocity data comes fro m iTuratto et al.l (|l99a ) ; iJha et al. 
(|l999l ): iGarnavich etsll J^OoJ) and Ixaubenberger et al.l (|2008l l. 
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Fig. 1 3. — Cartoon SNe la distribution. The observed subgroup fraction data comes from 

(|2nid v 



